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We demonstrate that a Hermitian matrix model can be derived from level 
truncated open string field theory with Chan-Paton factors. The Hermitian 
matrix is coupled with a scalar and U{N) vectors which are responsible 
for the D-brane at the tachyon vacuum. Effective potential for the scalar is 
evaluated both for finite and large N. Increase of potential height is observed 
in both cases. The large N matrix integral is identihed with a system of N 
TIL branes and a ghost FZZT brane. 


1 Introduction 

Two different candidates for the nonperturbative formulation of string theory have been 
known — string field theory (SET) and matrix model. It is commonly believed that 
they are just different descriptions of the underlying theory. Therefore, it is important 
to investigate the relationship between these formulations. However, a few examples 
have been known to reduce open SET (OSFT) [T] to certain matrix models. First is 
OSFT for topological A or B model [2] which reduces to Chern-Simons matrix model [3] 
or ordinary Hermitian matrix model [1] respectively. Second is OSFT for (2,1) minimal 
string [H] which reduces to Kontsevich matrix model [B]. In addition, less direct examples 
for c = 1 m, c = 0 0 and critical [9| strings have been known. In those examples, 
each matrix model is obtained from different OSFT associated with particular boundary 
conformal field theory (BCFT). A systematic way to derive different matrix models from 
OSFT in fixed background has not yet known. Finding such method is important to 
study the background independence of OSFT. 

In this paper, we present an example for such method. Our idea is simple: instead of 
varying BCFT, we start with a general setup in critical string, and approximate the string 
field to its first few components. An example we will study is the level truncation in the 
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universal sector |10j of critical {D = 26) OSFT in which the string field is approximated 
at level n as 

^ =-00 + V'l +'!/'2 H- hV’nQ (1-1) 

Given this approximation, the OSFT action [1] immediately reduces to a matrix actiorj^ 

S = OmnTr [M^Mn] + femnpTr [MmMnMp] , (1.2) 

where is a Hermitian matrix. The approximation is known to work well in the first 
few level and improves quickly as the level increases m]- 

A possible reason that such method has not yet been examined is the lack of under¬ 
standing of Chan-Paton factors in OSFT. Although Chan-Paton factors can be intro¬ 
duced to OSFT consistently, their origin has not yet been explained. Recently, Frier 
and Maccaferri proposed a new construction of classical solutions of OSFT in terms of 
the regularized boundary conditions changing (BCC) projectors [16]. Their construc¬ 
tion covers a wide range of backgrounds including multiple D-branes, which is the main 
interest of present paper. Following their work, Kishimoto, Masuda, Takahashi and 
Takemoto (KMTT) demonstrated that N Chan-Paton factors naturally arise from the 
decomposition of string field in terms of the regularized BCC projectors m- String held 
theory expanded around the multiple D-branes solution was interpreted as a system of 
-|- 1 D-branes. We will employ their formulation as our foundation. 

This paper is organized as follows. Section[^introduces KMTT formulation. In section 
we derive a one-matrix model in terms of the level truncation. In section]^ we derive 
the partition function of the model both in hnite and large N. Section summarizes 
our results and further discussions are given. 

2 KMTT formulation 

Let us briehy review the KMTT formulation |17j of OSFT. They begin with Erler- 
Maccaferri solution for multiple D-branes [16j 

N 

= (2.3) 

a=l 

where is the tachyon vacuum solution and and Sq are the regularized BCC 
projectors which obey 


SaLlfe = 5ab, Qr'^a = = 0. (2.4) 

Here Qt is the kinetic operator at the tachyon vacuum dehned by 

= QbT-F + 4'4't- (2.5) 

^Each component field carries Chan-Paton indices. 

^ As we will see later, extra vectors and scalars couple with the matrix. Here we omit them for 
simplicity. 
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The index a in Sq or is a Chan-Paton factor which labels D-branes. The solution (2.3) 
carries degrees of freedom required for multiple D-branes. We further expect that non- 
Abelian gauge gauge symmetry is realized in OSFT expanded around (2.3). KMTT m 
introduced the decomposition of the string field which realizes this: 


'I' — X + Xa^a + ^aXa + '^a4>ab^b: 


( 2 . 6 ) 


where x, Xa, Xa and Xab are component string fields. By expanding OSFT action around 


'I'o according to (2.6), they derived the matrix action 


S = Si + S2 + 53 , 




52 = - 


^XQtX + gX" 


53 = -^ 


2XaQT0Xa + XaXXa + Xa 4 >abXb 


(2.7) 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


wher e Qto is the kinetic operator defined by Qto'^ = Qb '^ -|- 'I't'!'. The trace in 
(2.8) runs over indices of and identical indices in (2.10) are to be summed oveiQ 
KMTT [T7] claimed that (2.7) describes -|- 1 D-branes rather than N D-branes. In 
addition to the N unstable D-branes, there is a “D-brane” at the tachyon vacuum 
described by y. According to their interpretation, (j)ab connects two unstable D-branes a 
and b while Xa connects D-brane a with the D-brane at the tachyon vacuum, y represents 
fluctuation on the D-brane at the tachyon vacuum. 


A remarkable feature of the action (2.7) is the presence of the “vector” sector described 


by Xa and Xa, which has not been found in literature. The action S 3 is quadratic for 
the vectors therefore can be integrated out in the path integral. With assuming suitable 
gauge fixing, we perform path integral for Xa and Xa and obtain a determinant factor 


det((5ro + </> + x) 


-1 


( 2 . 11 ) 


in the partition function. In the following sections, we will evaluate this determinant by 
truncating the string field rather than imposing conventional gauge condition such as 
linear gauge [T8] . 


3 Truncation to matrices 

Although the string fields presented in previous section carry Chan-Paton factors, they 
also depend on infinitely many labels which distinguish each state in a BCFT. In order to 
make our analysis tractable, we introduce an approximation in which dynamical variables 
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We will ignore a constant shift in the action hereafter since it is not relevant for remaining discussions. 
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reduce to matrices. Let us consider the string field in the universal sectoijj truncated 
at level n, where the level is defined by eigenvalue of the kinetic operator of OSFT. 
Denoting the level k base element of the string held V’fcj we truncate each component in 
(2.6) as 


= = Xa = ^d^VA:, X = (3.12) 

k=l k=l k=l k=l 

where is a Hermite matrix, is a complex vector, is its complex conjugate, 
and is a real number. Note that these component helds do not depend on any other 
variables. As a concrete example of such truncation, we choose the expansion examined 
in Ref. [TH] for dressed Bq gauge [50] . In this case, the level k base is given by 




= cK^Bc 


1 

1 + K' 


(3.13) 


where K,B,c are the elements of the KBc algebra [5T]. We would like to study level 1 
truncation in which string helds are given by 


</>a6 = 

(3.14) 

Xa = 

(3.15) 

Xa = 

(3.16) 

X = 

(3.17) 


A matrix action obtained by this truncation is rather complex since it includes both 
level 0 and 1 helds. However, it is possible to reduce degrees of freedom further. It is 
know that level 1 helds cannot have cubic terms since |19j 



(3.18) 


Thus, level 1 helds are quadratic in the action and can be integrated out using equations 
of motion. With the help of the explicit values of products of the basis ■00 and V’l 
evaluated in m, a nontrivial solution of equations of motion for and 

turns out to be 


= 


r{i) 


(3.19) 


By this assignment, four string helds in ( |3.12 ) become proportional to Erler-Schnabl 
solution [20] : 

^t = c{1 + K)Bc-^. (3.20) 

i + A 


^The open string field in universal sector is build from the Virasoro generator L„ and conformal ghosts 
b„ and c„ on the SL{2,R) vacuum |10 |. 
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Thus, the four field^can be written as 

4>ab = Mab'^T-, Xa = Xa = Ca^T, X = i^T- 


(3.21) 


Although only few fields are included in the approximation, we expect that it well 
captures the essence of the dynamics of D-branes, as is the case in the conventional level 
truncation analysis. An evidence that supports our expectation is that the truncation 
interpolates between two analytic solutions: one is the Erler-Maccaferri’s multiple D- 
branes denoted as {Mat, ^a, ^a, t) = (0, 0, 0, 0) while the other is the perturbative vacuum, 
also denoted as {5ab, 0,0,-1). 

Let us derive the truncated action. The requirement that T-p reproduces the correct 
value of the D-brane tension is represented by 


"^tQb^t =-/ ^tQto^t = 0 , / ^tQt^t = H — 


1T^ 


TT^ 




vr^ 


Applying these to (|2.7|), we find 

1 


s = -T 


9 


/2 


Tr ( —-M^ + —M^j + + ^ait^ab + Mab)Cb 


where we have rescaled the open string coupling as 

1 3 


g/2 g2jj-2 


(3.22) 

(3.23) 

(3.24) 


For later convenience, we shift M to M + 1 and omit the prime in g'. Then, ignoring 
constant shift, we obtain an action 


S = -^[Ti W{M) + W{t) + ^(1 + t + M)e] , 


(3.25) 


where W{x) = In this way, the truncation (3.12) reduces OSFT action 

to a cubic matrix action coupled with a scalar and complex vector^ The partition 


function for (3.25) can be obtained by employing the standard technique of matrix 
modeQ As readily found in (3.25), the action is invariant under U{N) transformation 
M —)■ UMU\ ^ —>■ U^, ^ —>■ ^1/1, where U denotes an U{N) matrix. This is a residual 
gauge symmetry of KMTT action [2Tj. This U{N) symmetry can be fixed by diagonal¬ 
izing Mab to its eigenvalues Xa with inserting Vandermonde determinant {Xa — A;,)^ in 
the partition function. After performing Gaussian integral for we obtain a partition 
function 


N 


z = 


dXadt ( 


.-V 


(3.26) 


a=l 


®This truncation makes sense if we replace \1/t with another representation of the tachyon vacnum 
solution, since we do not require its explicit form in following analysis. 

® A matrix model with U(N) vector was proposed in |22| in a different setting. 

^ For example, see [23j. 
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where 


N 


N 


V = 


^(^) + ^j + log |i + Ad + 1| — log(AQ, — Aft)^. (3.27) 


a=l 


a=l 


a<b 


This can describes a system of iV +1 particles moving in the potential W. An eigenvalue 


Aq feels repulsive forces from other eigenvalues through the last term of (3.27). Also, it 


is attracted towards —t — 1 due to the second term in (3.27). 


It is interesting to compare our result with the matrix formulation of c < 1 noncritical 
string theories. It has been recognized that Hermitian one-matrix model serves nonper- 
turbative definition of (2, 2k + 1) minimal string theory. As an example, let us consider 
the action studied in Ref. 


(3.28) 


5 = -^ [TrlT(M)+V'(M-z)V’] , 


where -0 and il) are fermionic vectors. They can be integrated out so that insert a factor 


det(M — z) 


(3.29) 


in the matrix integral. In double scaling limit, this determinant is identified with a FZZT 
brane |26l ], and the fermionic vectors are identified with fermionic strings between a 
FZZT brane and a stuck of N ZZ branes [28]. In contrast, our action (3.25) contains 


bosonic vectors ^ and ^ rather than fermionic ones. Integration with respect to them 
yields a factor 

det(Mt-Fl)"^ (3.30) 

in the partition function. Comparison between ( |3.29 ) and (3.30) naturally identifies the 
determinant (3.30) as a ghost FZZT brane [2^ which cancels the effect of a FZZT brane 
1^ Unfortunately, corresponding observable in minimal string theory is not yet identified. 


4 Effective potential 

In this section, we will evaluate the effective potential for t in terms of saddle point 
equations for eigenvalues. In large N limit of t’Hooft expansion, saddle point configu¬ 
rations are leading contributions to the matrix integral. Even for finite N, saddle point 
configurations also offer a good approximation to the matrix integral when g is small. 


In either case, saddle point equations are obtained from the variation of (3.27): 

N 

W'(\ ^ J_ - _9 

^ o “T ''\a ~T~ -L 

N 


b^a 


9 


= 0, (a = l...A^), 


^ a=l 


1 


i + Aa -l-1 


= 0 . 


(4.31) 


(4.32) 


®The ghost D-brane has been proposed as an object which cancels the effects of a D-brane m ■ 
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4.1 N = 1 


Let us study the dynamics of the model at finite N. We begin with = 1, which 
corresponds to a system of an unstable D-brane and another D-brane at the tachyon 


vacuum. Denoting Xa as A, the potential (3.27) is given by 


D = 4 (Wit) + W{X)) + log |t + A + 1|, 

r 

(4.33) 

and saddle point equations are 



(4.34) 


(4.35) 


By combining above two equations, we obtain a coupling independent equation 

t(t + 1) = A(A + 1), (4.36) 


which has two roots X = t and A = —t — 1. The latter is not appropriate since it hits 
the singularity (f + A + 1)“^ in the partition function. Therefore, we choose t = X as our 
solution. Given this choice, (4.34) and (4.35) reduce to single equation 


1 / 2 ^ 1 


= 0 . 


This can be rewritten to a cubic equation, 

t[t + l)(2t + 1) + = 0 


(4.37) 


(4.38) 


whose roots correspond to saddle points. It is easily understood that while there are 
three roots for small g, two of them disappear beyond critical value of g. Let us discuss 
further details as follows. Small g expansion of these roots is 

-g^-3g^ + 0{g^), -1 - g^ + 3g^ + 0{g^), -^ + 2g^ + 0{g^). (4.39) 

It is useful to show positions of these roots in a plot of the potential W. Fig. and are 
such plots for different values of g. Fig. [T] shows that saddle points for small g; saddle 
points are placed within the region (—1,0) with equal intervals. As g become larger, 
two larger roots become closer and annihilate beyond the critical value 

= 2-1/23-3/4 _ 

Fig.§ shows a placement of a root after the annihilation. 

Analysis made here is merely a classical approximation and not full quantum treat¬ 
ment. Fortunately, we can integrate out one variable in the potential without relying on 
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Figure 1: Solutions at s' = 0.10 


W 



Figure 2: Solution at ^ = 1.0. 
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the saddle point approximation since present example is enough simple. Let us change 
the variables as 


t “t" A 


Then, the partition function can be written as 


t — X 
2 


(4.41) 


Z = 


dudv- 
du- 


1 


2 u + l 
1 


|2u + 1| 



{-2W{u) - {2u + l)v‘^) 
W{u^ 


(4.42) 

(4.43) 


where we performed Gaussian integral for v in the second line. Thus, we obtain one 
dimensional effective potential 


14fr(n) 


-^W{u) + - log \2u + 1|. 
5 2 


(4.44) 


We note that the saddle point approximation corresponds to setting v = 
since we our choice of the saddle point, X = t corresponds to it. Therefore, 


0 in (4.42), 


2 1 

14addle(^t) = -^Wiu) + -log|2u+ 1|. 

2 


(4.45) 


By comparing (4.44) and (4.45), we hnd that they only differ in the coefficient of the 


logarithmic term which is negligible for small g. Thus, the saddle point approximation 
works well for small coupling. 

Finally, let us compare the shape of the potential for different values of g. Fig. [^is a 
plot of (4.44) for small and large g. It is observed that the stable vacuum around tt = 0 
disappears for large g due to dominance of the logarithmic term. This phenomenon 
corresponds to the annihilation of saddle points which has already been observed in 
the saddle point analysis. Such dependence of the effective potential on the coupling is 
consistent with the fact that the D-brane describes the system at small coupling. 


4.2 N>2 


Next we proceed to > 2, where the equations of motions are given by (4.31) and (4.32) 


We again perform saddle point approximation by finding roots of these equations. The 
roots can be obtained numerically for each value of g. Numerical solutions for g = 0.01 
and up to iV = 3 are shown in tables Hi and We pick real roots only, and each 
solution is specified as (t, Ai,..., Aat). Due to the symmetry which exchanges As, we 
only need to specify single configuration for each value of t among N permutations of 
A. We also evaluate the value of effective potential V{X,t) normalized by the value of 
D-brane tension 1/(6^^). We observe an interesting pattern in these results. Each root 
is located at one of a ‘site’ which is obtained by dividing the region — 1 < A < 0 into 
-|- 1 intervals, i.e., 


m 


^m. — 


Tv + r 


(m = 0,1,... A^ -|- 1). 


(4.46) 
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Figure 3: Effective potential for u, where solid line is a plot for g = Q.l while dashed line 
is for g = 1.0. 


t 

A 

V 

-1.0001 

-1.0001 

2.0000 

-0.4799 

-0.4799 

0.9947 

-0.0001 

-0.0001 

0.0000 


Table 1: Numerical solutions for N = 1 sX g = 0.01. 


t 

Ai 

A 2 

A permutations 

V 

-0.6664 

-0.9998 

-0.3330 

2 

1.9950 

-0.3331 

-0.6664 

0.0098 

2 

0.9950 

-0.0002 

0.0098, 

-0.0102 

2 

0.0052 


Table 2: Numerical solutions for N = 2 aX g = 0.01. 


t 

Ai 

A 2 

A 3 

A permutations 

V 

-0.4997 

-0.9996 

-0.4999 

-0.0003 

6 

1.9955 

-0.2498 

-0.7496 

-0.0100 

-0.0099 

6 

1.0005 

-0.0003 

-0.0176 

-0.0169 

-0.0002 

6 

0.0155 


Table 3: Numerical solutions for N = ‘i aX g = 0.01. 
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We also observe that the minimum value of t is given by —2/{N + 1) while values of As 
can reach —1. We also observed that the sum of all eigenvalues are close to its value of 
potential, i.e., 

N 

t + (4.47) 

a=l 

Finally, it is also interesting to see that the maximum value of the potential is close to 
2 for all N. 

We would like to close this subsection with a summary of our result: 

• The stable vacuum is lost for large coupling. This indicates breakdown of D-brane 
description. The critical value of the coupling can be determined by saddle point 
equations. 

• At small g, the local maximum of the effective potential is about twice higher than 
that of original one. 


4.3 Large N 


Let us evaluate the partition function in large N limit following with the standard 
method of matrix model |23) . We begin with the partition function 


, na<b(^a - WiK)) 

ria + Aa + 1| 


and introduce the t’Hooft coupling 


Then, saddle point equations read 


^ = g^N. 


N . 2\ 

' rt 


1 


= 0 , 


Aa + i + 1 
Ai , 2\ 1 2 


= 0 . 


(4.48) 


(4.49) 


(4.50) 


(4.51) 


In large N limit, eigenvalues are described by a continuous distribution p(X) and sum¬ 
mation for eigenvalues is replaced with integration 


1 

N 


E 


d\p{X). 


(4.52) 


Then, saddle point equations are replaced with 


1 


{t +1^) + 


dX 


pW 

t + A -|-1 


= 0 , 


(4.53) 
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-/"S- 


(4.54) 


At leading order, the second term in (4.54) becomes negligible. Thus we obtain saddle 
point equations 

pW 


— f dX 

^ J 

-(A + A^)-2 [ dX 
M J 


A + t + 1 

. d(A0 
A-A' 


= 0, 

(4.55) 

= 0. 

(4.56) 


The latter equation (4.56) solves the planar limit of a cubic matrix model whose solution 


(4.57) 


can be found elsewhere [30]. It is convenient to introduce the resolvent 

pW 


uj{z) = dX 


X — z 


The equation (4.56) can be replaced with an equation for uj{z). Once uj{z) is obtained, 

(4.58) 


(4.55) can be solved by Ending solutions of 


1 


+ +w(-t - 1) = 0. 

The “one-cut” solution for the resolvent in large N limit is known to be m 
^ - a){z - b) + ^) } ’ 


(4.59) 


where a and b are endpoints of the branch cut, which define a support for the eigenvalue 
density p{X). Requiring u}{z) ~ z~^ at infinity, we obtain equations 


3(a + 6)^ + 4(a + 6) — 4o6 = 0, 

(a + b)^ + (a + 6)^ — 4a6(a + 6 + 1) = 16//. 


(4.60) 

(4.61) 


which solve a and b as functions of //. These equations are conveniently rewritten in 
terms of parameters a = a + b and a = a — b as 


16// = -cr((T + l)(cj + 2), 
2 {a + lf + a = 2. 


(4.62) 

(4.63) 


The latter equation restricts a inside —2 < cr < 0. Further, (4.62) tells us that there are 

(4.64) 


no real roots of (4.62) within —2 < cr < 1. Therefore a is constrained within 

- 1 < C7 < 0. 


Let us choose a brunch which starts from <7 = 0. From (4.62), the maximum value of // 
reads 


cTc = ^(-3 + \/3)~-0.42, 


1 


Pc = 


24\/3 


0.024 


(4.65) 
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Given these ingredients, equation (4.55) can be rewritten into 


- I -(t + + -^t‘^ + {a + 2)t+ -0-2 + 2cr + 1 | ~ 0. 


Integrating this equation yields an effective potential for t, 

V{t) = -(Ul(t) + V2{t)) 


where 


viit) = -W{t), 


V 2 {t) = ^ y f z‘^ + {a + 2)z + ^ 0-2 + 2(7 + 1. 


(4.66) 


(4.67) 


(4.68) 


(4.69) 


A plot of effective potential (4.67) is shown in Fig. . The local minimum at t = 0 and 


V 



Figure 4: Effective potential at large N. The value of the potential is normalized by 
l/(6/r). The original potential W{t) is shown together. 

the local maximum around t = —1 are observed for any value of /r. This is consistent 
with the t’Hooft limit, in which g is small so that the D-brane description of the system 
holds. On the other hand, the height of the local maximum is always higher than that of 
the original potential W. The height is maximum for ^ = 0; in this case, the potential is 
twice higher than W. While the height decreases as g increases, it remains higher than 
W even at a maximum value of /r = /ic- The increase of the height from the original one 
can be understood from the particle description of eigenvalues. Eigenvalues filled in the 
bottom of W{t) pull the tachyon t by attractive force. It makes harder for t to climb 
the potential wall, thus increases the height of the potential height. Such dependence 
of the effective potential on /r is quite different from that of the probe eigenvalue model 
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investigated in m where a plateau along the eigenvalue distribution is observed. The 
existence of the plateau is explained by the fact that a probe eigenvalue cannot be 
distinguished from others. Our model has no plateau since t can be distinguished from 
other eigenvalues. 


5 Conclusions and discussions 


In this paper, we proposed a systematic method to derive matrix models from level 
truncated OSFT. Obtained matrix model contains U{N) vectors and a scalar in addition 
to Hermite matrix. We have evaluated the effective potential of the scalar both for finite 
and large N. Increase of the potential height was observed at small coupling. In section 
we have interpreted our model as a system of ZZ branes and a ghost FZZT brane. 
We would like to discuss further issues to be explored. First, we would like to present 
an alternative but rather heuristic interpretation of our result. Let us go back to the 
inverse of the determinant : 

det{QTo + (p + x)~^ ■ (5.70) 

The basic idea is that this quantity can be regarded as a propagator with (j) and x 
insertions. Recall that this factor is obtained by an integrating out Xa and Xa which 
connect a D-brane with the tachyon vacuum where the world-sheet boundary disappears. 
Therefore, it is natural to think that (5.70) amounts to a disk amplitude with single 


world-sheet boundary of a D-brane. Schematically, such disk amplitude can be written 


Tr / —e~dQTo+<t>+x) _ 


(5.71) 


As is well known, small t limit of such amplitude corresponds to closed string propaga¬ 
tion [32]. Therefore, the inverse determinant (5.70) encodes gravitational force between 
N D-branes and the tachyon vacuurrj^ 

Second issue is abut the level truncation. We have seen that the approximated OSFT 
action at first few levels yields a one-matrix model which can be interpreted as c < 1 
noncritical string theory [33| . It is interesting to improve the approximation by including 
higher level fields to obtain multi-matrix models. We speculate that that the improved 
matrix action continues to be dual to some closed string theory on more nontrivial 
background. Finally, we will recover original OSFT with infinitely many matrices in 
infinite level limit. We also expect that this matrix model describes a critical string 
theory in nontrivial background through AdS/CFT like duality |34j . Thus, the level 
truncated OSFT offers a way to describe closed strings in approximated geometry. 

Last issue is the matrix description of OSFT based on the left-right splitting of open 
strings which have been examined in past []. Our model based on the KMTT decompo¬ 
sition looks quite differently from these models. However, as mentioned in the previous 
paragraph, our model will recover full OSFT in infinite level. Thus the left-right split¬ 
ting models and our model both describe same OSFT. We expect that all ingredients of 


®This interpretation is consistent with the attractive force between A and t observed in section 
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KMTT, including Chan-Paton factors and regularized BCC projectors, are embedded 
in the left-right type matrix model in quite nontrivial manner. 

Together with recent development which deals with different backgrounds as classical 
solutions m, our result presents further evidence for SFT as a formulation of nonper- 
turbative string theory. We hope that further developments in this direction will shed 
light on the landscape of string theory. 
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